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ABSTRACT 

The Mount Wilson Ca HK survey revealed magnetic activity variations in a large sample of solar-type stars 
with timescales ranging from 2.5 to 25 years. This broad range of cycle periods is thought to reflect differences 
in the rotational properties and the depths of the surface convection zones for stars with various masses and 
ages. In 2007 we initiated a long-term monitoring campaign of Ca II H and K emission for a sample of 
57 southern solar-type stars to measure their magnetic activity cycles and their rotational properties when 
possible. We report the discovery of a 1.6-year magnetic activity cycle in the exoplanet host star t Horologii, 
and we obtain an estimate of the rotation period that is consistent with Hyades membership. This is the shortest 
activity cycle so far measured for a solar-type star, and may be related to the short-timescale magnetic variations 
recently identified in the Sun and HD 49933 from helio- and asteroseismic measurements. Future asteroseismic 
observations of i Hor can be compared to those obtained near the magnetic minimum in 2006 to search for 
cycle-induced shifts in the oscillation frequencies. If such short activity cycles are common in F stars, then 
NASA's Kepler mission should observe their effects in many of its long-term asteroseismic targets. 
Subject headings: stars: activity — stars: chromospheres — stars: individual(HD 17051, HR 810) — surveys 



1. BACKGROUND 

Astronomers have been making telescopic observations of 
sunspots since the time of Galileo, gradually building a histor- 
ical record showing a periodic rise and fall in the number of 
sunspots every 1 1 years. We now know that sunspots are re- 
gions with a sufficiently strong magnetic field to alter the local 
thermal structure, so this 1 1 -year sunspot cycle actually traces 
a variation in surface magnetism. Attempts to understand this 
behavior theoretically often invoke a combination of differ- 
ential rotation, convection, and meridional flow to m odulate 
the gl obal field through a magnetic dynamo (e.g., see Rempel 
2006). Although we cannot observe spots on other solar-type 
stars directly these areas of concentrated magnetic field pro- 
duce, among other signatures, strong emission in the Ca II H 
(396.8 nm) and K (393.4 nm) spectral lines. The intensity 
of the emission scales with the amount of non-thermal heat- 
ing in the chromosphere, making these lines a useful spectro- 
scopic proxy for the strength of, an d fractional area covered 
by, magnetic fields (Leighton 1959J). IWilsonl (I 1978b was the 
first to demonstrate that many solar-type stars exhibit long- 
term cyclic variations in their Ca II H and K (hereafter Ca HK) 
emission, analogous to the solar variations. 

Significant progress in dynamo modeling emerged af- 
ter helioseismology provided meaningful co nstraints on 
the Sun's interior stru cture and dynamics dBrown et alj 
Il989t iSchou et ail ll998). Variations in the mean strength 
of the solar magnetic field lead to significant shifts 
(~0.5 /x Hz) in the frequencies of even the lowest-degre e 
p-modes (iLibbrecht & Woodardlfl990l ISalabertet all 12004). 
Space-based astero seismolog y missions, such as MOST 
(Walk er et al.l 120031) CoR oT dBaglin et al.1 12006b . and Ke- 
pler dBorucki et all 120101) . as well as ground-based net- 
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works like the Stellar Observations Network Group (SONG; 
iGrundahl et ail 120081) . are now allowing additional tests 
of dynamo models using other solar- type stars (e.g., see 
Chaplin e t alj|2007t IMetcalfe et aill2007l) . 

The F8V star i Horologii (i Hor = HD 17051 = HR810, 
V=5.4, B-V=0.57) hosts a non -transiting 2 Mj exoplanet with 
an orb ital period of 311 days dKurster et al.ll2000t iNaef et al l 
12001b . Although it is currently situated in the southern hemi- 
sphere, kinematic considerations have led to the suggestion 
that it could be an eva porated member of the Hyades clus- 
ter (M ontes et al.H200ll) . Asteroseismic observations support 
this conclusion, since the acoustic oscillation frequencies of 
the star are best reproduced with models that have the same 
metallicity, helium abundance, and stellar age as other Hyades 
members (IVauclair et al. 2008|). 

We report the discovery of a 1 .6-year magnetic activity cy- 
cle in l Hor from synoptic Ca HK measurements obtained 
with the Small and Moderate Aperture Research Telescope 
System (SMARTS) 1.5-m telescope at Cerro Tololo Inter- 
american Observatory (CTIO) since 2008. We provide an 
overview of the survey methodology and analysis methods 
in §|2] and we present the stellar activity measurements and 
other derived properties in |3] We conclude with a discus- 
sion of the broader implications of this discovery for stellar 
dynamo modeling and future observations in ^4] 



2. OBSERVATIONS & DATA REDUCTION 



The chromospheric activity survey of iHenrv et al.l ([1996) 
contained a total of 1016 observations of 815 individual stars 
with visual magnitudes between 0.0 and about 9.0, which 
were observed using the RC Spec instrument on the CTIO 
1.5-m telescope. Several sub-samples were defined, including 
the "Best & Brightest" (B) and "Nearby" (N) samples, which 
together contain 92 individual stars with visual magnitudes 
between 0.0 and 7.9, and B-V colors that are approximately 
solar. In August 2007, we began a long-term Ca HK monitor- 
ing program for the 57 sta rs in the combined (B +N) sample 
that are brighter than V= 6 (IMetcalfe et al.ll2009b . the limiting 
magnitude of future ground-based asteroseismic observations 
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by SONG. All of the most promising southern asteroseismic 
targets are included in this B+N sub-sample. 

Since January 2008, we have obtained 74 low-resolution 
(R ~ 2500) spectra of 1 Hor covering 37 distinct epochs with 
the upgraded RC Spec instrument on the SMARTS 1 .5-m tele- 
scope. Using standard IRAF0 routines, the 60 s integrations 
were subjected to the usual bias and fiat field corrections and 
the spectra were extracted and wavelength calibrated using 
a reference He-Ar spectrum o btained immediately before the 
stellar exposures. Following iDuncan et alj d!991h . the cali- 
brated spectra were then integrated in 1.09 A triangular pass- 
bands centered on the cores of the Ca H and K lines and 
compared to 20 A continuum passbands from the wings of 
the lines to generate a CTIO chromospheric activity index, 
Sctio- These values were converted to Mount Wilson in- 
dices (Smwo) using data for 26 targets that were observed con- 
temporaneously with the Solar-Stellar Spectrograph at Low- 
ell Observatory (J. Hall, private communication), and adopt- 
ing a quadratic function for the correlation (cf. He nry et alj 
1996). The details of our observations of l Hor are listed in 
Table Q] Note that the uncertainties shown in Table Q] repre- 
sent only the internal errors. The systematic uncertainty from 
the conversion between the CTIO and MWO activity indices 
is cr sys - +0.007 (MWO - CTIO). 

5 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy (AURA) under cooperative agreement with the National Science Foun- 
dation. 



In addition to the single-epoch observation from 
IHenrv et all (fl996h . 5 M wo = 0.225 ± 0.005 on 1992.9479, 
there are several other Ca HK measurements of 1 Hor in 
the literature that w e can use to probe a ctivity variations on 
longer timescales. Jenk ins et alj (2006) publish ed a recal- 
ibration of measurements originally made by iTinnev et al. 
(2002) on 2001.5918 with S M wo = 0.249 ±0.002. iGrav et al. 
(2006) obtained a spectrum on 2002.9538 with a revised 
Smwo = 0-226 ± 0.01 (R. G ray, private communication). Fi- 
nally, ISchroder et all (120091) measured S M wo = 0.246 ± 0.03 
from a spectrum obtained on 2003.9387. Additional Ca HK 
measurements of l Hor have appeared in the literature, but 
without sufficient detail to determine the precise epoch of the 
observation and the S index on the Mount Wilson scale. 

3. RESULTS 

Our time-series measurements of the Mount Wilson S index 
for l Hor are plotted in FigureQ] We fit a sinusoid to these data 
and found a cycle period of f cyc =1.6 years around a mean 
value of (S) = 0.242. The variations are not expected to be 
strictly sinusoidal, so the fitted amplitude of A cyc = 0.024 does 
not capture the full range of observed values from S ~ 0.21- 
0.28, nearly 30 % of the mean activity level (the solar range 
is S ~ 0.17-0.20; Baliunas et a O 19951) . The parameter values 
of the fit do not change significantl y when incl uding the few 
archival data points, and all but the IGrav et ail measurement 
agree with the extrapolated sinusoid at the la level. 

We were initially intrigued that the activity maximum at 
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FIG . 1 . — Chromospheric activity measurem ents of the F8V star i Hor from 
the southern HK survey (Metcalfe et al. 2009), showing a clear variation with 
a cycle period of 1.6 years, the shortest cycle measured for a sun-like star. 
Note that the error bars represent only the measurement errors, and do not 
include the systematic uncertainty <T sys ~ 0.007 (arrow). 



2009.3 coincided with an epoch of periastron for the eccen- 
tric exoplanet in the system (cf. iNaef et al.l l2001). Although 
activity i nduced by hot Jupiters has been seen in a few cases 
(e.g., see Walker et al. 2008|, we would not expect it in this 
case because even at periastron the star-planet separation is 
0.7 AU. The exoplanet orbital period of 3 1 1 days bears no ob- 
vious relation to the cycle period of 584 days, so this appears 
to be simply a coincidence. 

Aste roseismic observations of l Hor by IVauclair et alj 
(2008) in 2006.9 coincided with the magnetic minimum one 
cycle prior to the beginning of our data set. As first observed 
in the Sun more than two decades ago (Libbrecht & Woodard 
1990), the global oscillation frequencies are shifted sig- 
nificantly from solar minimum to maximum. The ampli- 
tude of these frequency shifts was predicte d to be larger 
(~1 uHz) for star s hotter than the Sun (see iMetcalfe et al.l 
120071 IKaroffetal.i r2009). as was re c ently confirmed in the 
F5V star HP 49933 by iGarcfa et ail (|2010). Future astero- 
seismic observations near the magnetic maxima that coincide 
with the observing season for i Hor in late 2010 or 2013 have 
the best chance of detecting these frequency shifts. 

The scatter around the sinusoidal variation in Figure [TJ is 
caused in part by rotational modulation of individual active 
regions. Based on prior single-epoch measurements of its 
chromospheric activity le vel, the rotation peri od of l Hor was 
estimated to be 7.9 days (Saar & Osten 1997). After normal- 
izing our measurements by the 1.6 year sinusoid, we passed 
the residuals throug h a Lomb-Scargle periodogram dLombl 
Il976t IScargld U982) to search for the signature of rotation. 
The results are shown in Figure [2] with the highest peak at 
f rot = 8.5 ±0.1 days and a smaller peak near 7.9 days. We 
verified that a single sinusoid with a period of 8 days, sampled 
in the same way as our data, produces a clear and significant 
peak in the periodogram. Of course, the rotational modula- 
tion of individual active regions will not be coherent over the 
span of our data set, and spots at different latitudes may have 
different periods. These effects will tend to reduce the sig- 
nificance of individ ual peaks in the pe riodogram. However, 
a recent analysis by Bois se et all (120 10) of the radial velocity 
measurements from IVauclair et al.l (12008b found evidence of 
rotation in t Hor with periods in the range 7.9 to 8.4 days, con- 




8 10 12 

Period (days) 

FIG. 2. — The Lomb-Scargle periodogram of our Ca HK measurements 
after removing the 1 .6 year sinusoid, suggesting a rotation period that is con- 
sistent with Hyades membership for this star. 



sistent with the dominant peaks in Figure [2] from our Ca HK 
measurements. The mean rot ation period of sol ar- type stars 
in the Hyades is ~8.4 days dRadick et alj|1995l). so both of 
the se results support the conclusions of Montes et alj (1200 lb 
and IVauclair et al.l ((2008) that i Hor is an evaporated member 
of the cluster. 

If we combine the asteroseismic radius of i Hor, R = 
1.18 Rq (S. Vauclair, private communication) and the rota- 
tion period above with th e measured v sin i = 6 .47 ±0.5 km s" 1 
from Butler et al. (2006), we can obtain an estimate of the in- 
clination angle, i ~ 60°. If the rotation axis is perpendicular 
to the orbital plane of the exoplanefj then the absolute mass 
of the planet is only about 15% larger than the minimum mass 
obtained from the radial velocity orbital solution. 

4. CONCLUSIONS & DISCUSSION 

The immediate question that arises from the observation of 
such a short activity cycle is whether we can understand it in 
the context of dynamo models that are frequently invoked to 
explain the 1 1 -year solar cycle. It is generally difficult to ex- 
trapolate solar dynamo models to other stars because it is not 
well understood how the underlying parametrizations change 
for different stellar properties. However, in the case of t Hor 
we can take advantage of a lucky coincidence — an F8 star 
with a rotation period near 8 days ha s a Coriolis number 2£It c 
that is very similar to that of the Sun dKiiker & Riidigerl2005l) . 
Since the Coriolis number characterizes the influence of rota- 
tion on convection, it determines to a large extent the overall 
direction and profile of the turbulent Reynolds stresses (RS) 
that drive differential rotation and meridional flow. The am- 
plitudes are expected to be different due to the higher lumi- 
nosity and resulting convective velocities (RS ~ v^), but it still 
allows us to make some simple estimates of how various dy- 
namo scenarios considered for the Sun are expected to scale 
to i Hor. 

From a standard stellar evolution model of l Hor 
with the global pro perties inferred from asteroseismology 
(IVauclair et alj |2008). we deduce that the convective velocity 
v c is 2.8 times higher than the solar value while the relevant 

6 Note that recent measurements of the Rossiter-McLaughlin effect for 
some transiting hot Jupiter systems suggest a substantial spin-orbit misalign- 
ment (e.g., see Schlaufman 2010). 
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timescale r c is 3.2 times shorter (due to a smaller pressure 
scale height). Using simple scaling relations these quantities 
allow us to estimate turbulent diffusivities T) t , v, that are 2.4 
times solar. The meridional flow speed v m is expected to be 
about 2.6 times solar (from angular momentum transport bal- 
ance, v 2 c ~ v m O) and the differential rotation 2-3.3 times solar 
(considering thermal balance £!A£! ~ AT and angular mo- 
mentum transport balance v 2 c ~ i/ t AfT). Additional scaling 
factors of order unity can arise from the different depth and 
aspect ratio of the convection zone. Nevertheless these esti- 
mates agree very well with the differential rot ation model for 
an F8 star presented in Kiiker & Riidigei] (l2005h . 

The above properties would lead us to predict a cycle pe- 
riod of 4-5 years for an advection dominated flux-transport 
dynamo (P cyc ~ v" 1 ), or 3.5-5 years for a classical a-Vl dy- 
namo (P cyc ~ [Aila] -0 ' 5 , assuming a ~ v c ). These estimates 
are both more than a factor of two longer than the observed 
cycle, but for the flux transport dynamo we assumed that the 
underlying conveyor belt would extend to similar latitudes as 
observed in the Sun. Due to the larger aspect ratio of the con- 
vection zone in an F8 star as well as the larger magnetic diffu- 
sivity, it is conceivable that the flux transport is cut short lead- 
ing to a significant reduction in the dynamo period. Currently, 
B. P. Brown et al. (in preparation) are exploring dynamo ac- 
tion in solar-type stars rotating at 3-5 times the solar rate us- 
ing global 3D anelastic MHD models. For a sufficiently large 
level of turbulence, cyclic behavior is found with periods of a 
few years or less. These studies demonstrate that computing 
a detailed 3D model of i Hor should already be possible. 

Prior to the discovery of this 1.6-year magnetic cycle in 
i Hor, the shortest measured cycle periods were 2.52 years 
(HD 76151) and 2.60 years (HP 1 90406) from the Mount Wil- 
son survey (Baliuna s et aljfl995l) . Both of these appeared to 
be sec ondary cycles super imposed on a much longer primary 
cycle. Bohm- Vi tensd d2007l) has suggested, based on the sam- 
ple of stars with well characterize d rotation and cycle peri- 
ods in Saar & Brandenburg ( 1999), that the "active" and "in- 
active" branches in the Prot-^cyc diagram may be caused by 
two distinct dynamos that are driven in different regions of 
the star. Specifically, iBohm- Vitensd suggests that the active 
branch may represent a dynamo operating in the near surface 
shear layer, while the inactive branch is driven by the shear 
layer at the base of the convection zone. In this scenario, stars 
that exhibit cycle periods on both branches must have the two 
types of dynamos operating simultaneously. If the 1.6-year 
cycle in i Hor is on the inactive branch, then the active branch 
cycle period is expected to be near 6 years. Although the pre- 
vious Ca HK measurements are sparse, we see no evidence 
of a secular trend in Smwo on longer timescales. Continued 
observations by our program should yield stronger constraints 
on possible slow variation s. 

Recently, iGarcfa et al.l ( 12010b detected the signature of a 
short magnetic activity cycle in the F5V star HD 49933 us- 
ing asteroseismic measurements from the CoR oT satellite 
( Appo urchaux et al]l2008t iBenomar et al J 120091) . Just as in 
the Sun, where the global oscillation modes shift to higher 
frequencies and lower amplitudes towards the maximum of 
the 11-year solar cycle, HD 49933 appeared to pass through 
a magnetic minimum during 137 days of continuous observa- 
tions. An additional 60 days of data from an earlier epoch 
could not place strong constraints on the cycle period, but 



suggested a value between 120 days and about 1 year. If con- 
firmed by ground-based monitoring of the Ca HK lines, this 
would place HD 49933 (with P mt = 3 .4 d) in the same category 
of magnetic cycle observed in i Hor. 

The asteroseismic signature of a short magnetic cycle has 
also been detected in the Sun itself. iFletcher et alJ (2010) 
analyzed the low-degree solar oscillation frequencies from 
the BiSON and GOLF experiments and found evidence of a 
quasi-biennial (2 year) signal in both data sets after removing 
the dominant 11 -year period. Unlike the 11 -year signal, the 
amplitude of the 2-year variation appeared to be largely inde- 
pendent of frequency, leading the authors to suggest that the 
secondary cycle must be operating independently. However, 
the amplitude of the 2-year signal was uniformly larger during 
the maximum of the 11 -year cycle, suggesting that buoyant 
magnetic flux might be rising from the base of the convection 
zone and pumping up a near-surface dynamo with the 2-year 
period. Active branch stars with a cycle period of 1 1 years 
are expected to show a secondary cycle period on the inac- 
tive branch around 2 years. However, this normally occurs in 
stars rotating at twice the solar rate, and the identification of 
the two dynamos would then be reversed — with the 11 -year 
dynamo operating in the near-surface shear layer, while the 
2-year dynamo is driven at the base of the convection zone. 

If short magnetic activity cycles are common, NASA's Ke- 
pler mission should detect them in the asteroseismic measure- 
ments of many additional stars. In principle such measure- 
ments can provide unique constraints on the underlying phys- 
ical mechanism, and Kepler will also yield measurements of 
some of the key dynamo ingredients. Even without the short 
cadence data for asteroseismology, the high precision time- 
series photometry from Kepler is sufficient to characterize the 
surface di fferential rotation through detailed spot modeling 
(e.g., see iBasri et al.l l20Tot) . For the brighter asteroseismic 
targets where the individual oscillation frequencies are de- 
tectable, the time series will be long enough to resolve ro- 
tational splitting of the modes into multiplets for stars with 
rotation rates betw een about 2 and 10 times the solar rate 
dBallot et al.ll2.008b . Measurements of the rotational splitting 
as a function of radial order can indirectly probe radial differ- 
ential rotation, since the various modes sample slightly differ- 
ent (but overlapping) regions of the star. For the very best and 
brightest asteroseismic targets, Kepler will obtain a frequency 
precision sufficient to measure the depth of the surface con- 
vection zone from the oscillatory si gnal in the so-calle d sec- 
ond frequency differences (e.g., see IVerner et al.ll2~006h . The 
Kepler mission is expected to document all of these properties 
in at least a few dozen solar-type stars, gradually leading to a 
broader context for our understanding of the solar dynamo. 
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